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Abstract. In this paper we describe in detail the reduction, preparation and reliability of the photometric cata- 
logues which comprise the CFH12K- VIRMOS deep field. This region, consisting of four contiguous pointings of 
the CFH12K wide-field mosaic camera, covers a total area of 1.2 deg 2 . The survey reaches a limiting magnitude 
of Bab ~ 26.5, Vab ~ 26.2, Rab ~ 25.9 and Iab ~ 25.0 (corresponding to the point at which our recovery rate 
for simulated point-like sources sources falls below 50%). In total the survey contains 90,729 extended sources in 
the magnitude range 18.0 < Iab < 24.0. We demonstrate our catalogues are free from systematic biases and are 
complete and reliable down these limits. By comparing our galaxy number counts to previous wide-field CCD 
surveys, we estimate that the upper limit on bin-to-bin systematic photometric errors for the I— limited sample 
is ~ 10% in this magnitude range. We estimate that 68% of the catalogues sources have absolute per co-ordinate 
astrometric uncertainties less than |Aa| ~ 0.38" and A<5 ~ 0.32". Our internal (filter-to-filter) per co-ordinate 
astrometric uncertainties are |Aa| ~ 0.08" and |A<5| ~ 0.08". We quantify the completeness of our survey in the 
joint space defined by object total magnitude and peak surface brightness. We also demonstrate that no signifi- 
cant positional incompleteness effects are present in our catalogues to Iab < 24.0. Finally, we present numerous 
comparisons between our catalogues and published literature data: galaxy and star counts, galaxy and stellar 
colours, and the clustering of both point-like and extended populations. In all cases our measurements are in 
excellent agreement with literature data to Iab < 24.0. This combination of depth and areal coverage makes 
this multi-colour catalogue a solid foundation to select galaxies for follow-up spectroscopy with VIMOS on the 
ESO-VLT and a unique database to study the formation and evolution of the faint galaxy population to z ~ 1 
and beyond. 

Key words, observations: galaxies - general - astronomical data bases: surveys: cosmology: large- scale structure of 
Universe 



1. Introduction 

Deep multi-colour observations remain one of the sim- 
plest and most useful tools at our disposal to investigate 
the evolution and properties of the faint galaxy popu- 



lation. The first digital studies of the distant Universe 


llMetcalfe et all 


1991; 


Lilly et all 


ll99lllTvson & Seitzerl 


1988 


pushed catalogue limiting 


magnitudes beyond the 
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traditional Bj = 24 boundary of previous photographic 
works. They provided the first glimpse of the Universe at 
intermediate redshifts and t entative upper limits on th e 
numbers of z > 3 galaxies (jGuhathakurta et all fl990i. 
Following these studies, t housand-object spectroscopic 
survey s in the mid- nineties (Ellis et all ll99fiHFillv et~aTL 
1995.3), combined with reliable photometry provided by 



charge-coupled devices (CCDs) mounted on 4m tele- 
scopes, sketched out for the first time the broad outlines 
of galaxy and structure evolution f or normal L* galaxies 
in the redshift ra nge < z < 1 l|Le Fevre et all Il996t 
iLillv et allll995blh 

At higher redshift s, the Hubble Deep Fields project 
(Wil liams et al.. 1996) has demonstrated just how much 
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can be accomplished beyond the spectroscopic limit of 
even 10m telescopes, provided photometric observations 
are sufficiently accurate. And in the last five years, search- 
ing for the passage of the Lyman break through broad- 
band filters has proved to be an ext remely efficient way 
to locate objects at z > 3 and above llSteidel et allll996t 
iMadau et aD. Il996t ISteidel fc HamiltorJ.ll993|) . These re- 
cent studies have provided a unique insight into the forma- 
tion of large-scale structures and the relationship between 
mass and light in the regime where theoretical predictions 
depend sensitively on th e underlying model assumptions 
ijAdelberger et alj . Il998|) . Meanwhile , at lower redshifts, 
the 2dF, 2MASS and Sloan surveys jcolless etafl. l200ll 
IJarrett et all 1200(1 lYoTk" et all l2000j) are now providing 
us with a detailed picture of local Universe and a solid 
reference point for studies of the distant Universe. 

To progress beyond these works at z ~ 1 demands 
imaging catalogues at least an order of magnitude larger 
than previously available from surveys based on single 
CCD detectors. One simple reason is that many ob- 
jects of interest are rare; surface densities of Lyman- 
break galaxies, for instance, are around ~ 1 arcmin - 
at Iab < 24.0; amassing large samples of such galax- 
ies requires coverage of a correspondingly large area 
of sky. In the recent years, large format multiple-CCD 
cameras such as the Universi ty of Hawaii's 8K camera 
(UH8K; iLuppino et aD l|l994|) 'l. containing many sepa- 
rate CCD detectors and with a field of view of around 
~ 0.25 deg 2 have become available and several surveys 
are now underway or have been recently completed us- 
ing these cameras. For example, the Canada-France Deep 
Fields Survey (CFDF) covered three surv ey fields of 
the C anada-France Redshift Survey (CFRS; iLillv et all 
1995a) with the UH8K camera, prov iding UBVI pho- 
tome try for around 50, 000 objects ijMcCracken et all 
With even larger, more efficient instruments such as 
the Canada-France-Hawaii Telescope's 12K mosaic cam- 
era (CFH12K) it has now become feasible to conduct 
surveys subtending degree scales on the sky and reach- 
ing mean redshifts of z ~ 1. Furthermore, the advent 
of very high throughp ut spectrographs such as VIMOS 
( Le Fevre et all l200f)|) demands a new generation of wide- 
area deep digital survey catalogues in which systematic 
errors are rigorously controlled. It is also important to 
develop robust and reliable tools which can effectively 
deal with the even larger volumes of data which will be 
produced by the forthcoming 500-night Canada-France- 
Hawaii- Telescope Legacy Survey (CFHTLS) project which 
will u se the new 1 deg 2 fiel d-of-view MEGACAM instru- 
ment l|Boulade et all l2000|) . 

In this paper we describe BVRI photometric obser- 
vations which have been carried out in the VIRMOS- 
VLT deep field, a 1.2 deg 2 high galactic latitude field 
of low galactic extinction which will be the centrepiece 
of the VIRMOS-VLT deep spectroscopic survey (VVDS; 
iLe Fevre et "all 1200 1|) . Several tens of thousands of galaxy 
spectra will be acquired here to a limiting magnitude 
of Iab < 24.0. Observations described in this paper 



were made using the CFH12K mosaic camera: the broad 
spectral response and sensitivity of 12K, compared to 
previous-generation cameras such as the UH8K, make it 
an ideal instru ment for this kind of s urvey (the DEEP2 
redshift survey I Davis fc FaberL Il998h are also using the 
CFH12K camera for survey pre-selection) . In addition to 
the core BVRI dataset described here, near-ultraviolet, 
near-infrared, X-ray and 1.4 GHz radio observations have 
also been carried out in this field, and will be presented 
in forthcoming papers. A companion paper (Le Fevre et 
al., paper I) describes the overall strategy of the survey. 
Finally, a third paper will describe shallower optical data. 
This comprises three 4 deg 2 fields (the fields at 14hrs, 
22hrs and lOhrs) and a 2 deg 2 extension to the field de- 
scribed here. In each of these fields spectra will be ac- 
quired from a sample of galaxies magnitude limited at 
Iab = 22.5. 

This paper is organised as follows. In Section [2] we 
describe observations which were carried out; Section |21 
describes how this data was reduced and includes de- 
tails of how both astrometric and photometric solutions 
were computed to produce the final output images. In 
Section 13.51 we explain how object catalogues were pro- 
duced from these images. In Section 0] we present a de- 
tailed quality assessment of the catalogues, where com- 
parisons are made to existing deep catalogues. Our con- 
clusions and summary are presented in Section |SJ 

2. Observations 

This paper describes the F02 deep field of the 
VIRMOS-VLT survey. This field is centred at 02' 1 26"W 
-04°30'00" (J2000) which corresponds to a galactic lat- 
itude and longitude of (172.0, —58.0) respectively. All 
observations described here were carried out using the 
CFH12K camera over a series of three observing runs from 
November 1999 to October 2000. The CFH12K is a 12,288 
by 8192 pixel mosaic camera comprising twelve 2048 x 4096 
15 micron pixel thinned backside illuminated MIT Lincoln 
Laboratori es CCID20 devices. It h as a total field of view 

IH[2000). The average gain over 



of 42' x 28' ijCuillandre 
all CCDs is 1.6e~/ADU. The pixel scale of the detector is 
approximately 0.205"/pixel. The CFHT wide-field correc- 
tor produces a radial distor tion amounting to ~ 15 pixels 
at the camera field of view ( Cuil landre et al. L ll996h which 
necessitates resampling our images before stacking. In to- 
tal, the F02 deep field comprises four pointings of CFH12K 
in four filters covering a total area of 1.2 deg 2 . 

Cosmetic quality of the detectors is generally good, 
with less than 200 bad columns in total, most of which 
are concentrated on a single detector. In addition, CFHT 
provides pixel masks which allow these bad pixels to be 
easily removed. The average gap between the charge cou- 
pled devices (CCDs) in the east west-direction is 7.8", 
whereas the north-south gap is 6.8". To fill these gaps 
in the final stacked mosaic, exposures were arranged in a 
dither pattern of radius 15". Figure Q shows the layout of 
all /—band observations made on the field. Each rectan- 
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Fig. 1. An outline of all CCD frames contained in the I— 
band stack. Shown in this Figure are four separate point- 
ings of the CFH12K camera. At each pointing there are 
7-8 individual exposures comprising 12 CCDs in each. At 
the scale of this plot, it is not possible to see the displace- 
ment caused by the dither pattern offsets. 
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Fig. 2. Seeing distribution for CCD images used in the 
F02 mosaic for each of the four filters. The dashed line 
shows the median seeing; note that in all cases this is less 
than 1". All images with a seeing higher than 1.2" were 
discarded and are not shown in this Figure. 



gular box corresponds to the outline of one CCD image. 
At the scale of this plot it is not possible to see the dither 
pattern displacements. 

Exposure times per image in each filter were typi- 
cally 1800s for B, 1080s for V, 1200s for R and 720s 
for /. The total exposure times in each of the four fil- 
ters at each of the four pointings are listed in Table ^ 
In Figure [21 we show the distribution of full- width half- 
maximum (FWHM) seeing values for all the input images. 
All data with FWHM values greater than 1.2" were dis- 
carded. 



3. Data reductions 

Reducing observations taken by mosaic camera detectors 
present unique challenges in terms of the quantity of data 
which must be processed. For the project described in 
this paper, several thousand separate CCD images were 
processed totalling several hundred gigabytes of data; at 
least double this quantity of disk space is needed to com- 
bine these images into a single mosaic. The very large 
volume of data involved meant that all processing was 
carried out at TERAPIX, a data reduction centre at the 
Institut d'Astrophysique de Paris designed specially to 
process large quantities of wide-field mosaic camera data. 
One of the principal objectives of the TERAPIX facility is 
to carry out high-level processing (image co-addition and 
catalogue extraction) for the CFHTLS survey. The soft- 
ware tools described here were developed in the context of 



this project 1 , and are described in detail in the following 
Sections of this paper. These tools were designed to be 
as general as possible and will work on mosaic data from 
other telescopes and cameras. They are available at the 
TERAPIX web site 2 . 

Our data reduction technique involves first de-trending 
our raw data, then computing astrometric and photomet- 
ric solutions for all input CCD frames. These solutions 
are then used to re-sample and stack all CCDs for each 
filter into a single contiguous mosaic, from which cata- 
logues are then extracted. We now describe each of these 
steps in turn. 

3.1. Prereductions 

All data were pre-reduced using the Fits Large Image 
Processing Software (FLIPS) package 3 . Pre-reductions 
followed the normal steps of overscan, bias and dark 
subtraction. Dark current for the CFH12K is negligible, 
around 3 e~pix _1 hr . Following this, all images are then 
divided by a twilight flat. For data taken in R and I fil- 
ters, fringe removal is accomplished by subtracting a 2a 
clipped median super-flat. These superflats were typically 
constructed from around 60 science images. At this stage 
all images are scaled so they have the same sky back- 
ground. After pre-reductions, the gradient of the back- 

1 http:/ /terapix. iap.fr/doc/doc. html 

2 http://terapix.iap.fr/soft 

3 http://www.cfht.hawaii.edu/ jcc/Flips/flips.html 
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Field R.A. (2000) Dec. (2000) Band Exposure time Median seeing Completeness limit 



(Pf, P2,P3,P4) (hours) (arcsec) (AB mags) 



0226-0430 02:26:00 -04:30:00 B 6.0, 5.5, 5.5, 5.0 0.9 26.5 

V 4.0, 7.0, 7.0, 4.0 0.8 26.2 

R 3.8, 3.1, 3.7, 3.3 0.8 25.9 

I 3.8, 3.2, 3.0, 2.8 0.8 25.0 



Table 1. Details of the images used. We list the total integration time in each of the four pointings and the completeness 
limit of each filter, based on simulations described fully in Section l4~Tl 



ground variation is measured across the CCD frames. This 
is always less than 2% and is typically ~ 0.5%. 

3.2. Computing the astrometric solution 

Astrometry and relative photometry were computed using 
the "astrometrix" and "photometrix" tools, part of the 
"WIFIX" package for the reduction of wide-field images 4 , 
developed in the context of the TERAPIX project. In the 
following sections we describe the algorithms underlying 
these tasks and how they were applied to our data. An 
important aspect of these procedures is that they are able 
to operate in an unattended mode; given the large volume 
of data involved, manual interactions must be kept to a 
minimum. 

In the following discussion, the approximate area de- 
fined by one of our CFH12K frames is described as a 
'pointing'; in the F02 field there are four pointings in to- 
tal. To reach a specified depth in a given pointing, several 
images are required, and these exposures are offset from 
each other in a "dither sequence" which contains between 
five and ten exposures. 

To co-add all these separate CCD images and produce 
a seamless output image requires a detailed knowledge 
of the astrometric solution of each CCD frame. However, 
computing an astrometric solution constrained only by 
an external astrometric catalogue will not produce a suf- 
ficiently accurate solution. This is because the intrinsic 
accuracy of even the best external catalogues have rms 
uncertainties of around ~ 0.2" per co-ordinate. This level 
of accuracy (combined with low object sky densities) make 
astrometric solutions derived from these catalogues alone 
unsuitable for co-addition of frames. In our survey we wish 
to produce a single contiguous mosaic comprising four sep- 
arate pointings in four separate bands. Overlapping im- 
ages from different filters and adjacent pointings must be 
registered to sub-pixel accuracy (or an r.m.s. smaller than 
~ 0.1". For these reasons, we have developed a proce- 
dure in which the astrometric solution is constrained si- 
multaneously by the external astrometric catalogue and 
an internal catalogue generated from matched objects in 
overlapping CCDs, both from objects in the same pointing 
and in adjacent pointings. Since our I— band data gener- 
ally have the best seeing, and the VIRMOS spectroscopic 
sample will be I— selected , we compute an astrometric 



http:/ /www. na.astro.it/~radovich/wifix. htm 



solution from the I— band data first. We then use this 
stacked image to derive solutions for all other bands. 

The astrometric catalogue we us ed is the United States 
Naval Observatory (USNO)-A2.0 l(Moned . Il998h which 
provides the position of ~ 0.5 x 10 8 sources. After the re- 
moval of saturated and extended objects between 30 and 
50 sources per CFH12K CCD frame are left (the distri- 
butio n of astrometric sources is not homogeneous on the 
skv). lDeutschl lll999h . by carrying out comparisons with 
a sample of ~ 300 radio sources, estimate a radial posi- 
tional accuracy of the USNO-A2.0 of ~ 0.25" (defined as 
the radius enclosing 68% of the objects). 

We begin our astrometric procedure by extract- 
ing input cat alogues from the (J-ba nd data) using 
SExtractor iBertin fc Arnoutsl (|l99f)() on each input 
(flat-fielded and bias subtracted) image. Saturated sources 
were removed. A cross-correlation procedure was used to 
compute the initial offsets (generated by the dither se- 
quence) between these catalogues and the astrometric ref- 
erence catalogue. 

The astrometric solution describes the transformation 
between pixel projection plane co-ordinates (x, y) and 
spherical co-ordinates on the sky £, 77, namely, 



i=Q j=Q 
p i 

i=0 j=0 

This is equivalent to the TAN projection in 
ICalabretta fc Greisenl l)2002|) but without the (x 2 + y 2 )°' 5 
terms. 

This solution is computed as follows: firstly, for each 
observing run we select an exposure for which the relative 
offsets between CCDs, rotation angles and scaling factors 
are computed. These values are then applied to all other 
exposures. In this way it is possible to use all the avail- 
able data to compute the offset to the USNO. In order 
to further minimise the possibility of errors, we separate 
the exposures into groups where fields overlap by around 
50% and then compute the offsets between exposures in 
the same group. The offset to the USNO must be the same 
for a given group and we take the median of these values. 

Next, for each CCD we match image sources with the 
astrometric catalogue and recompute the linear part of 
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Fig. 3. Positions of USNO sources compared to point-like 
sources with 17.5 < Iab < 22.0 in the final 1.2 deg 2 
stacked image. The inner and outer circles (at 0.5" and 
1.0") enclose 68% and 90% of all objects respectively. The 
dotted lines cross each other at the position of the cen- 
troid, which is at (Aa, AS) = (-0.03", 0.00"). 

the world coordinate system (WCS). We then compute a 
first-guess astrometric solution (Eq^with p = 2). At this 
stage, a database of overlapping sources from all CCDs 
and pointings ("master" catalogues) is also constructed. 

In the final step we re-compute the astrometric solu- 
tion for all CCDs, but this time including points both 
from the astrometric catalogues and from the overlap- 
ping sources catalogue. The equation is solved using as 
weights the positional uncertainties on each source: we 
assume 0.3" for the USNO sources, whereas for the over- 
lapping sources we use the positional uncertainties given 
by SExtractor. The system of equations is solved by 
iteration. Three iterations with p = 3 produces RMS per- 
coordinate uncertainties of ~ 0.3" for the USNO sources 
and ~ 0.05" for overlapping sources. Unfortunately, since 
the offsets between exposures in the same dithering se- 
quence are small (< 20 — 30"), the astrometric solution is 
poorly constrained in some regions as overlapping sources 
are actually extracted from the same part of the CCD and 
are affected by the same distortion (discussed at length 
later in this section). We therefore decided to reduce the 
weighting of overlapping sources from exposures which 
come from the same dither sequence. For these objects, 
weights are divided by the number of frames in that dither. 

How do the positions of objects in the final I— stack 
(produced after resampling and co-adding of all the in- 
put CCDs using this astrometric solution) compare to 
our input USNO-A2 catalogue? In Figure |3] we plot the 
difference between the positions of non-saturated, point- 
like sources in the I— stack and their counterparts in the 




Fig. 4. Positions of USNO sources compared to point-like 
objects in the / stack as a function of right ascension and 
declination, displaying the uncertainty per co-ordinate. 
The inner dashed lines enclose 68% of all sources and are 
at |Aa| - 0.38" and |A<5| - 0.32"; the outer dot-dashed 
lines enclosing 90% of all sources are at |Aa| ~ 0.80" and 
|A<5| ~ 0.70". 



USNO-A2. The inner and outer circles enclose 68% and 
90% of all the objects respectively, and have radii of 0.5" 
and 1.0". Figure 01 plots these residuals as a function of 
right ascension and declination. We note that our use of 
overlapping sources ensures that the accuracy of our in- 
ternal astrometric solution is much higher, which is nec- 
essary if we are to successfully produce coadded images, 
as we will discuss now. 

To compute the astrometric solution for the BVR 
data, we use as our reference catalogue a list of non- 
saturated point-like sources extracted from the I— band 
stack. The procedure described above is then repeated us- 
ing this catalogue. The surface density of these sources 
is much higher than that of the USNO-A2, and this ap- 
proach allows images with different filters to be registered 
to a much higher accuracy. In Figure El we show the po- 
sitional uncertainties between objects extracted in the B 
and / 1.2 deg 2 stacks. The inner circle and outer circles 
plotted at 0.11" and 0.26" enclose 68% and 90% of all ob- 
jects. These figures indicate we have reached our required 
level of sub-pixel accuracy. In Figure these positional 
uncertainties are plotted as a function of right ascension 
and declination. 

Examining Figure^l however, a potential problem be- 
comes apparent. It is clear that there are regions in the 
stacked images where the rms astrometric uncertainty is 
much higher the average, although the area affected by 
these problems is small, less than 1% of the total area. 
Two effects are responsible for these difficulties. The first 
is that the surface density of astrometric reference sources 
on the sky is not homogeneous: a situation can easily arise 
for example, in which there are few reference sources on 
one half of a CCD. Secondly, the dither offsets we use are 
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Fig. 5. Positions of unsaturated point-like sources in the I 
stack image compared with objects extracted from the B 
stacked image. For clarity only one point in ten is shown. 
The circles at 0.11" and 0.26" enclose 68% and 90% of 
all sources respectively. The dotted lines cross each other 
at the position of the centroid, which is at (Aa, A5) = 
(-0.01", -0.01"). 
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Fig. 6. Residuals between I and B stacks as a function 
of right ascension and declination. The inner dashed lines 
enclose 68% of all sources and are at |Aa| ~ 0.08" and 
|A5| ~ 0.08"; the outer dot-dashed lines enclosing 90% of 
all sources are at |Aa| ~ 0.18" and |A<5| - 0.18". 



too small to adequately sample the astrometric solution 
of the detector, and in many cases overlapping sources 
come from the same regions of the same CCDs. The net 
result is that in certain areas (mostly at the edges of CCD 
frames) the astrometric residuals are large. However, these 
regions are easily identified by comparing stacked images 
from each filter, and these areas are masked. The amount 



of area lost from these problems is small, less than a few 
percent of the total. 

3.3. Photometric calibration 

The photometric calibration is carried out in two steps. 
In the first, we determine an absolute calibration based 
on observations of standard star fields. In the second, we 
combine the absolute photometric calibration with a rela- 
tive calibration computed using catalogues of overlapping 
sources to derive scaling factors which will be applied dur- 
ing stacking. We now describe both these steps in turn. 

The absolute photometric calibration was performe d 
by comparisons to a set of Landolt fields llLandoltlll99l . 
We were able to identify at least one set of photometric 
images from our observations, based on local sky condi- 
tions and stability of the object fluxes. Our observations 
are scaled to these photometric images using the proce- 
dure outlined below. Note that in this work, we do not 
attempt to calibrate each individual CCD frame indepen- 
dently, but instead calibrate the scaled, flat-fielded images 
(i.e., the entire mosaic). In other words, we assume that 
each CCD has the same colour equation. Since our ob- 
servations were made, the CFH12K photometric perfor- 
mance has been extensively characterised as part of the 
CFHT queue observing program Elixir, which has deter- 
mined zero-points and colour equations based on obser- 
vations of several hundred standard star frames 5 . The as- 
sumption that each CCD has the same photometric equa- 
tion produces an overall photometric residual of less than 
4%. These investigations also indicate that the CFH12K 
BVR filters have negligible colour terms (less than 0.03 in 
(B — V) and (V — R) with respect to the Johnson Kron- 
Cousins system). Only in / was a significant colour term 
found (quoting from the ELIXIR web pages): 

he = ic f hi2k + 26.185 - 0.04(k - 1) + 0.107(i? - I), (2) 

where k represents the airmass, R, I are Johnson mag- 
nitudes and Ike is the magnitude in Johnson Kron-Cousins 
system. 

These zero-points and colour equations are broadly 
consistent with those determined from our survey data, 
but are of higher accuracy. Note that we do not actually 
use these colour equations in this work, and present all 
our magnitudes in the CFH12K instrumental system, de- 
fined by the combination of CFH12K detectors, filters and 
telescope. We do however convert our inst rumental m ag- 
nitudes to the "AB" magnitude system, ((OkeL Il974h as 
follows: B AB = b cfhl2k - 0.094, V AB = v cfhl2k - 0.007, 
Rab = r cf hi2k + 0.218 and I AB = i c fhi2k + 0.433, com- 
puted from the CFH12K filter curves. Our zero-points are 
corrected for galactic extinction of E(B — V) = 0.027 es- 
tim ated at the centre of the field provided by the maps 
of ISchlegel et al.1 i|l998h . As the galactic extinction is 
quite uniform over the field, we do not attempt to apply 
position-dependent magnitude corrections to our data. 

5 http:/ /www.cfht. hawaii.edu/Instruments/Elixir/filters. html 
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After computing the absolute photometric calibration, 
and identifying photometric exposures in our data, we 
next compute the relative scaling exposure-to-exposure. 
The zero-point of each image i may be written as the sum 
of two components: 



Z, = zP h 



(3) 



where the photometric zero point z ph is determined using 
standard stars and the relative zero point z\ incorporates 
the airmass correction and other changes in the atmo- 
spheric transparency. We assume that the relative zero 
point is the same for all the CCDs in a mosaic and that at 
least one image has been taken in photometric conditions, 
so that it may be used as reference. It is computed by min- 
imising the average difference in magnitude of overlapping 
so urces in different fram es . Our approach is similar to that 
of iKoranvi et al.l l)l998j) except that we do not derive our 
coefficients by direct solution of a matrix inversion but by 
an iterative minimisation process. 

The astrometric solution detailed in Section 13.21 pro- 
vides a database of overlapping sources which can be used 
to compute the photometric offsets. We first extract from 
these catalogues pairs of objects from overlapping frames 
with magnitudes mi,rrij and photometric uncertainties 
u>i,Wj. For each pair we compute 



Y^Wimi E' 



E 



(4) 



Our aim is to find the set of values z,; which minimises 

5m*j —< to >i +Zi - (< to >j +zj) = Srriij + (z t — Zj)(5) 

For each frame i with overlapping sources in other iV ov 
frames, the zero point is: 



Z t = z ph + 



E 

1 _r_ L ~'3 



N 



Zj — Srriij 



E" ph 

iVph : 



N 



(6) 



where we assume that iVph frames are known to be pho- 
tometric and that for them the average relative zero point 
is null. The effective airmass k, of the coadded image is: 



Ef ph ^ 



P h 



(7) 



all input images) using the astrometric solutions computed 
in Section EOl We use a tangent-plane projection, which 
is adequate given the relatively small size of our field. The 
flux scalings discussed in the previous Section are also ap- 
plied at this stage. An interpolated, projected weight-map 
is also computed which incorporates information derived 
from bad pixel maps. Both for the weight-map and image 
data, interpolation is carried out using a "Lanczos-3" in- 
terpolation kernel. This kernel corresponds to a sine func- 
tion multiplied by a windowing function. It provides an 
optimal balance between preserving the input signal and 
noise structure and limiting ringing artifacts on image dis- 
continuities (for instance, around bright stars and satellite 
trails). Additionally, at this stage a local sky background 
is computed at each pixel and subtracted, which ensures 
that the stacked image is free of large-scale background 
gradients. 

In the final stage, all resampled images are coadded 
(with each pixel weighted by the appropriate weight map) 
to produce the final output image. Pixels are coadded us- 
ing a median, which although sub-optimal for signal-to- 
noise properties, provides the best rejection of satellite 
trails and other cosmetic defects. In addition to the fi- 
nal coadded image, a total weight map is also produced 
containing information concerning how often each pixel 
has been observed. This weight-map is used during the 
detection and catalogue generation steps, outlined in the 
following sections. An important aspect of the resampling 
and stacking stage is that a single, large contiguous im- 
age is produced which contains all data, resulting in only 
one catalogue for each filter. No merging of catalogues is 
required. We set the size, pixel scale and orientation of 
the BVR data to match the one computed for the / data 
so that catalogues are easily extracted from these other 
images. 

3.5. Catalogue preparation 

To prepare merged BVRI catalogues we locate objects 
using a separate detection image and then perform pho- 
tometry in each bandpass at the positions defined by the 
detection image. Thi s image is construc ted using the x 2 
technique outlined in lSzalav et al.l l|l999l) . namely 



V N (5m* ) 2 

Eq. EI is iterated until — - — — converges or a maxi- 
mum number of iterations (50) is reached. 



(8) 



3.4. Image resampling, stacking and flux scaling 

Once the astrometric and photometric solutions have been 
computed for all the input CCDs, these images are com- 
bined to produce the final stack. This is carried out in 
a two-ste p process using a new image resampling tool, 
SWarp 6 (iBertin et al.ll2002t . 

In the first step, each input CCD is resampled and pro- 
jected onto a subsection of the output frame (which covers 

6 http:/ /terapix. iap.fr/soft/swarp 



where dj represents the background-subtracted pixel 
value in filter i, Gi the rms noise at this pixel and n is the 
number of filters (four in this case). 

In the absence of sources, the distribution of sky val- 
ues in this stacked image can be described as \ 2 distri- 
bution with n degrees of freedom. Detection and thresh- 
olding, therefore, have relatively simple interpretation in 
probabilistic terms. Our primary motivation in using this 
technique is to simplify the generation of multi-band cat- 
alogues and to reduce the numbers of spurious detections. 
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We produce this \ 2 image from the four stacked im- 
ages using SWarp which performs the local estimations 
of Oi and the sky background at each pixel (note that in 
version 1.21. of SWarp, used here, the output pixels in 
the generated x 2 images are actually S — \J (y/n)). For 
the x 2 technique to work reliably, seeing variations across 
input images must be small (less than ~ 20%) and more- 
over there must be the same number of input images at 
each pixel (a consequence of this is that detections may 
not be reliable in parts of the image where there is no cov- 
erage in all four bands). We have compared magnitudes 
of objects recovered with and without the use of the x 2 
detection image and we find no evidence that it introduces 
any magnitude bias. 

For objects to be included in our catalogues, they must 
contain at least three contiguous pixels with S > 2 in 
the x 2 image. This corresponds to a per-pixel detection 
threshold of 2.8a. With this threshold we detect in to- 
tal 460,447 sources. We chose this somewhat conservative 
threshold in order to minimise the numbers of false detec- 
tions. However, as we shall see in the following sections, 
our catalogues are essentially complete at Iab — 24. This 
threshold is close to the 'optimal' value where the proba- 
bility of misclassifying object pixels as s ky and sky pi xels 
as objects is minimised (see Figure 1. of lSzalav et al.h . 

In our catalogues, magnitudes are measured using 

iKronl l)l980(l -like elliptical aperture magnitudes computed 
using SExtractor (the magauto parameter). These 
magnitudes provide the most reliable measurement of the 
object's total flux (based on simulations carried out on 
our images), although they may be subject to blending 
effects if there are near neighbours. This is not a concern 
for us because our images in all bandpasses are far from 
the confusion limit and the numbers of objects affected 
are very small. We set a minimum Kron radius of 1.2", 
which means that for faint unresolved objects (where the 
Kron radius can be difficult to reliably determine), our 
"total" magnitudes will revert to simple aperture magni- 
tudes. Colours are always measured at the same location 
(defined by the detection image) on each stack. Because 
there are many large, extended objects on our images, we 
also use magauto magnitudes to measure object colours. 
Once again, tests show that these colours revert to "aper- 
ture" colours for the majority of faint, unresolved objects 
in our catalogues. Measurements show that for total mag- 
nitudes in the range 18.0 < Iab < 24.0 the difference 
between the median galaxy colour using total and aper- 
ture magnitudes is less than 0.05 magnitudes. 

We use the SExtractor f luxjradius parameter to 
identify point-like sources, using object profiles measured 
in the /—band images with centroids derived from the % 2 
image (of course, given that we have BVRI photometry 
for all objects, a complementary approach would be to 
identify stellar sources based on best-fitting spectral tem- 
plates; application of this technique to our data will be de- 
scribed in future papers). This parameter, denoted as ri, 
measures the radius which encloses a specified fraction of 
the object's total flux (in this case, 50%). Fainter than sat- 



uration (Iab ~ 18), point-like objects have a flux radius 
independent of magnitude. In our data, an added com- 
plication is that the seeing differs by around 20% across 
the mosaic, which means this locus is slightly dispersed. 
To account for this, we used an "adaptive" classification 
technique. In this procedure, we divide the field into many 
sub-areas and compute the mode of ri for objects within 
each box. Since point-like objects all have the same half- 
light radius this provides a robust measurement of the 
local seeing value. A search is then carried out for the bin 
which has a value of less than 20% of the modal bin; all 
objects with ri smaller than this are classified as stars. 

Figure|7|demonstrates how this method works in prac- 
tice. In the upper panel we show the half-light radii values 
as a function of Iab magnitude for all objects, normalised 
by the approximate mean value for the full field. In the 
lower panel, each object's half-light radii has been divided 
by the local measurement. The heavy dots in both panels 
indicate the point-like objects located using the classifier. 
The separation between resolved and unresolved sources 
is considerably improved, and distinguishing between the 
two seems feasible until Iab ~ 22, although in this work 
we adopt a conservative limit of Iab ~ 21.5. Beyond this 
limit no classification is attempted. 

Finally, the catalogue was visually inspected by over- 
plotting on the images all non-saturated objects with 
Iab < 24.0. Areas contaminated by residual fringing, 
satellite trails and scattered light were flagged by drawing 
polygonal masks around them. All large bright stars and 
their corresponding diffraction spikes were also masked. In 
Figure[H]we show the distribution of galaxies with 18.0 < 
Iab < 21.0. Masked regions, outlined by polygons, are 
also indicated. In the magnitude range 18.0 < Iab < 24.0 
we detect 90,951 extended objects not in the masked re- 
gions or at the borders of the image. In this area and 
magnitude range there are 2857 point-like sources. The 
total area of our survey, after masking is 1.18 deg 2 . 

4. Data quality assessment 

In these sections we describe and present a series of qual- 
ity assessment tests carried out on the BVRI catalogues 
prepared in the previous sections. 

4.1. Surface brightness selection effects and catalogue 
limiting magnitudes 

Our objective here is to determine the photometric com- 
pleteness of our survey in the plane defined by surface 
brightness and total magnitude. The importance of sur- 
face brightness selection effects in photome tric surveys 
has already been dis cussed by several authors l|Lillv et all 
ll995aHYoshulll99^ . In this work, an important objective 
of our simulations is to ensure that our data is fully com- 
plete and free from surface brightness selection effects at 
Iab = 24.0, the spectroscopic limit of the VLT-VIRMOS 
deep survey. 
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Fig. 7. Object half-light radius as a function of Iab mag- 
nitude. In the upper panel, we show the half-light radius in 
pixels (normalised by full-field mean of 2.5 pixels) over the 
full mosaic both for point-like (solid points) and extended 
sources (light points). In the lower panel each object's 
half-light radius has been normalised by the local mode, 
computed using the adaptive classifier (in both panels, for 
clarity, only a fraction of all objects are shown). 

To estimate how surface brightness selection effects 
could affect our catalogues a variety of approaches may 
be used. For example, we might try to produce simulated 
objects with a variety of light profiles, such as exponen- 
tial disks or r 1 / 4 -law bulges. These objects could then be 
dimmed based on some physically motivated picture of 
how galaxies evolve and then coadded back into our image 
frames and their object parameters re-measured. In what 
follows, however, we adopt a purely empirical approach. 
We start by extracting BVRI sub-areas of 100 bright, iso- 
lated galaxies from a 1024 x 1024 subsection of the field. 
Our objective is to quantify, at each combination of mag- 
nitude and total surface brightness, how many objects are 
recovered by our selection procedure. Some of the surface- 
brightness/magnitude combinations are not physical, but 
we wished to introduce as few assumptions as possible 
about the evolution of the faint galaxy population into our 
simulations. A consequence of this approach, of course, is 
that interpreting the overall completeness requires addi- 
tional information concerning where the faint galaxy pop- 
ulation resides in the surface-brightness/magnitude plane. 
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Fig. 8. Distribution of bright galaxies. The points show 
galaxies (extended sources) with 18.0 < Iab < 21.0. 
Masked regions are also indicated. 



To relocate our galaxies in the surface bright- 
ness/magnitude plane we reduce their total magnitude 
and surface brightnesses by dividing each object by some 
number to scale it to a fainter total magnitude. To reduce 
the surface brightness of the galaxies while keeping the 
same total magnitude, the galaxies were stretched using 
bilinear interpolation. After these alterations, the galaxies 
are added back into the subsections extracted from each 
of four bandpasses. Following this, objects are detected 
using the same procedures as in the real catalogues. We 
generate a \ 2 detection sub-image and search for objects 
using our detection threshold (corresponding to three con- 
nected pixels having S > 2). For each combination of Iab 
total magnitude and peak surface brightness (defined as 
the surface brightness of the brightest pixel in the recov- 
ered object) we measure the fraction of galaxies retrieved. 
To evaluate the numbers of spurious detections, all the 
images are multiplied by — 1. As before, the x 2 image is 
generated and the flux measurement is carried out using 
SExtractor in dual image mode on the negative im- 
ages. The number of initial detections will be the same 
since the chi-squared detection image generated from the 
negative images is the same as the chi-squared image gen- 
erated from the positive images. However, since the mea- 
surement image is now negative, positive fluxes will be 
measured only for the false detections. 

The results of these simulations are illustrated in 
Figure where they are plotted in the plane defined by 
total magnitude and peak surface brightness. Filled dots 
show the position of all objects in the field; the solid line 
shows the location of the stellar locus. For each combi- 
nation of total magnitude and peak surface brightness, 
the numbers indicate the percentage of objects which are 
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recovered. The contour lines indicate the 50%, 70% and 
90% completeness limits. The open circles show the spu- 
rious detections; we note that these spurious detections 
only become significant faintwards of Iab ~ 25. 

For point-like stellar objects at Iab = 24.0 our sur- 
vey is essentially complete, with approximately 90% of 
objects recovered at this combination of magnitude and 
surface brightness. At this magnitude our recovery frac- 
tion drops to 50% for objects with p eak surface bright- 
nesses of around ^25.5 mag arcsec 2 . iLillv et all 1 1995ajl 
have computed the tracks of typical galaxies in the surface- 
brightness magnitude plane, and their Figure 7 shows 
that most "normal" galaxies would be easily detected 
by our survey. Objects with extremely large half-light 
radii, such as the prototypical "Malin 1" low surface 
brightness galaxy could fall below our detection limits 
at higher redshifts. However, space-based observations 
suggests that the numbers of such extended, low-surface 
brig htness objects at the s e magnitudes is actu ally quite 
low l(Simard et all Il99flt iFasano et all Il998l) . Most of 
the objects we detect in the real data (filled points in 
Figure [5J fall in quite a narrow locus in the surface- 
brightness/magnitude plane. Moreover, if we compare our 
total galaxy number counts at Iab — 24.0 with an average 
of measurements in N orth and South Hubble Deep Fields 
( Me tcalfe et all 120011) (presented in Section 14. 3|) we see 
that we are in good agreement at Iab = 24 with these 
much deeper counts. Based on these considerations, we 
conclude that our survey is essentially free of selection ef- 
fects until at least Iab ~ 25. Finally, we note that our 
simulations also show that at this limit our survey is free 
from spurious detections (represented as the open circles 
in Figure |5J). 

We have also carried out a set of simulations based on 
adding progressively fainter simulated stellar sources to a 
sub-section of the image stacks. The magnitudes at which 
we recover 50% of the input sources is 26.5,26.2,25.9 and 
25.0 for BVRI filters respectively. 

4.2. Catalogue completeness as a function of position 

An important issue to address is how homogeneous the 
detection and measurement limits are over the entire mo- 
saic, and how stable the images' noise properties are as a 
function of location. To estimate this, we carried out a se- 
ries of simulations aimed at measuring the incompleteness 
as a function of position for the full I— band mosaic. In 
these simulations, we extracted ~ 200 sub-areas arranged 
in a grid over the mosaic. For each of these sub-areas, we 
added artificial stars at random positions using as a local 
seeing value the measurements derived in Section [3.51 four 
intention in this experiment is not to provide a realistic 
assessment of the absolute value of the incompleteness but 
rather its positional variation, which is why we adopt the 
simplifying assumption of stellar profiles). Detection and 
measurement was then carried out using SExtractor, 
using detections weighted by the local detection weight- 
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Fig. 9. Peak surface brightness as a function of Iab total 
magnitude. Detected galaxies are shown as the black dots, 
the open circles show the spurious detections. The stellar 
locus is represented by the solid line. Most galaxies lie 
nearby to this locus. The contour lines indicate the 50%, 
70% and 90% completeness limits. The numbers represent 
the fraction of objects recovered at each combination of 
total magnitude and peak surface brightness. 



map. For several slices in magnitude, we measured what 
fraction of stars we recover at each position. 

Our results are illustrated in Figure ^| where we plot 
the recovery fraction as a function of position for input ob- 
jects in the magnitude range 23.75 < Iab < 24.25. Over 
most of the field, the recovery fraction is > 90% (in good 
agreement with the result derived for point-like sources in 
the previous section). There are small regions where the 
completeness can be significantly lower; upon investiga- 
tion, these areas appear to be populated by bright stars 
which were not excluded from the simulation. From these 
tests we conclude that there is no significant variation in 
completeness over the mosaic at Iab = 24.0. 

4.3. Galaxy and star counts 

In this Section we compare the source counts of our cat- 
alogues of galaxies and stars (in reality, point-like and 
extended sources) to published compilations. The filled 
squares in Figure [21 shows -Z^-selected stellar num- 
ber counts compar ed to the predictions of the model of 
iRob in et al. I <l2002h . computed at the galactic latitude 
of our survey. In the magnitude range 18 < Iab < 
21.5 we measure the slope of the stellar number counts, 
d\og N/dm, as ~ 0.14. Our stellar counts and slope mea- 
surement are in in excellent agreement with model predic- 
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Fig. 10. Measured incompleteness as a function of posi- 
tion (decimal declination and right ascension) for the I 
band mosaic. The grey-scale contours indicate, at each po- 
sition, the fraction of stellar objects which were recovered. 
Input magnitudes were distributed uniformly in the range 
23.75 < Iab < 24.25. The thick lines show the location of 
the masked regions. 

tions (a detailed discussion of the properties of the stars 
in the survey will be presented in a future paper). 

Differential galaxy number counts per square degree 
per half magnitude for B, V, R, and I filters derived 
from our catalogue are represented as the filled circles in 
Figures El El El and El These counts are also listed 
in Table El We do not correct our counts for incomplete- 
ness, and show "raw" measurements. For our literature 
comparison, we include only measurements from recent 
digital surveys and omit earlier photographic works. We 
note how well the galaxy count measurements now agree 
well between the different surveys, in particular for the 
B— band measurements. It appears that the large scatter 
between different groups in previously presented compila- 
tions arose from uncertainties in photographic magnitude 
scales at the bright end and from the small fields of view of 
the previous generation of deep CCD surveys at the faint 
end. By contrast, in many large surveys today, system- 
atic errors arising from absolute zero-point uncertainties 
dominate in almost every magnitude bin. 

In general, our agreement with the literature measure- 
ments is excellent. In the V — band, our counts a re in 
good agreement with those of lCaban ac et al l ll200(tl. th e 
CFDF measurem ents and those of lArnoutsetalJ l|200l|) . 
The counts from IWilsonl l|2003h appear to be systemati- 
cally below all other measurements. 

We have also carried out a more detailed bin-by-bin 
comparison with the Iab galaxy counts determined in the 
CFDF survey, which has comparable total area (~ 1 deg 2 ) 
to the present work. The CFDF counts differ at most 
from the present work by ~ 10% per bin. Given the large 
numbers of galaxies in each bin, this error is almost al- 
ways larger than the combination of >/N Poisson count- 
ing errors and fluctuations induced by large scale struc- 



Filter 


Magnitude range 


dlogN/dm 


B 


20.0-24.0 


0.45 ±0.01 


B 


24.0-25.5 


0.32 ±0.04 


V 


20.0-24.0 


0.40 ± 0.03 


R 


20.0-24.0 


0.37 ±0.03 


I 


20.0-24.0 


0.34 ± 0.02 



Table 3. Best-fitting galaxy count slopes in BVRI bands 
for our catalogue, computed in the specified magnitude 
ranges. 
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Fig. 11. B^iB-selected galaxy counts (filled symbols), 
compared to a selection of recent results from the liter- 
ature. The dotted lines shows the best fitting slope com- 
puted in the magnitude ranges 20.0 < Bab < 24.0 and 
24.0 < Bab < 25.5. Error bars are smaller than our sym- 
bols and are not shown. 



ture (which is small for the fainter bins). Assuming this 
difference arises from uncertainties in the absolute zero- 
points of the CFDF and the present work, this would 
correspond to absolute zero-point uncertainty of ~ 0.05 
magnitudes (estimated simply by calculating what zero- 
point shift minimises the difference in counts between the 
two surveys). This is in agreement with what we have es- 
timated in Section 13.31 from observations of Landolt stan- 
dard star fields. 

We have also determined the slope a (dlogN/dm) of 
the galaxy counts, and these are reported in Tabled Our 
B and I slopes agree well with reported CFDF survey 
values of ~ 0.35 and ~ 0.47 respectively. Fainter than 
Bab ~ 24 we observe a slope flattening in the B— : for 
24 < B A b < 25.5, we find dlogN /dm — 0.3, similar to 
that found bv lMetcalfe et all l| 19951) . 
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B V R I 



AB Mag 


Ngal 


log N gal deg 2 


Ngal 


log Ngai deg 2 


Ngal 


log N ga i deg 2 


Ngal 


log N ga i de 


18.0-18.5 


37 


1.50 


137 


2.06 


303 


2.41 


332 


2.45 


18.5-19.0 


78 


1.82 


251 


2.33 


346 


2.47 


462 


2.59 


19.0-19.5 


119 


2.00 


335 


2.45 


508 


2.63 


811 


2.84 


19.5-20.0 


176 


2.17 


550 


2.67 


793 


2.83 


1156 


2.99 


20.0-20.5 


328 


2.44 


767 


2.81 


1157 


2.99 


1890 


3.20 


20.5-21.0 


508 


2.63 


1182 


3.00 


1732 


3.17 


2796 


3.37 


21.0-21.5 


882 


2.87 


1625 


3.14 


2476 


3.32 


4168 


3.55 


21.5-22.0 


1247 


3.02 


2565 


3.34 


3930 


3.52 


6759 


3.76 


22.0-22.5 


2312 


3.29 


4137 


3.54 


5922 


3.70 


9992 


3.93 


22.5-23.0 


3813 


3.51 


6831 


3.76 


9241 


3.89 


13676 


4.06 


23.0-23.5 


7014 


3.77 


11450 


3.99 


14006 


4.07 


19863 


4.23 


23.5-24.0 


12771 


4.03 


19609 


4.22 


22235 


4.28 


28828 


4.39 


24.0-24.5 


22130 


4.27 


30850 


4.42 


33632 


4.45 


39874 


4.53 


24.5-25.0 


34448 


4.47 


43593 


4.57 


46362 


4.59 


47079 


4.60 


25.0-25.5 


47255 


4.60 


54457 


4.66 


56686 


4.68 


43459 


4.57 



Table 2. Differential galaxy number counts. Counts are shown in half-magnitude bins selected in (BVRI)ab- Also 
shown are the logarithmic counts, normalised to the effective area of the survey field (1.18 deg 2 ). 
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Fig. 12. VAs-selected galaxy counts (filled symbols), com- 
pared to a selection of recent results from the literature. 
Error bars are smaller than our symbols and are not 
shown. 



4.4. Stellar colours 

In this Section we examine the relative and absolute pho- 
tometric calibration by investigating the colours of stars 
distributed over our field. 

The filled symbols in Figure fT~5l shows the (B — V)ab 
and (V — I)ab colours of point-like sources in our cat- 
alogue compared to synthetic colours generated by con- 



Fig. 13. it^s-selected galaxy counts (filled symbols), 
compared to a selection of recent results from the liter- 
ature. Error bars are smaller than our symbols and are 
not shown. 



volving the stellar spectral energy distributions of Pickles 
lll998|) with the CFH12K filter set (open symbols). From 
this plot it is apparent that in BVI filters our absolute 
calibration is accurate to within ~ 0.05 magnitudes. 

We also use measured colours of stars in our fields 
to investigate if there are systematic differences in stel- 
lar colours across the mosaic. The results of these tests 
are illustrated in Figures ED where we plot the median 
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Fig. 14. Iab — selected galaxy and star counts (filled cir- 
cles and squares respectively), compared to a selection of 
recent results from the literature (open symbols). For the 
galaxy counts, in all cases our (Poisson) error bars are 
smaller than the symbols and are not plotted. The dotted 
line shows the best fitting slope computed in the magni- 
tude range 20 < Iab < 24. The solid line shows the pre- 
dictions of the galactic model of Robin et al., computed 
at the galactic latitude of our survey. 




(v-i)„ 

Fig. 15. (B— V)ab vs (V-I)ab colours for (filled points) 
compared to synthetic colours computed from library of 
Pickles (open circles). 



Fig. 16. Stellar colours as a function of position for each 
of the four pointings. In each pointing we compute the me- 
dian (B — V)ab and (R-T)ab colours in bins of (V— I)ab 
(upper and lower panels respectively) . This is represented 
as the four different line styles. The regions are defined 
as follows: solid line: south-west, dotted line: north-east; 
dashed line: south-east, dot-dashed line: north-west. 

(B — V)ab, (V — I)ab and (R — I)ab colours of stars 
selected from each of the four pointings. It appears from 
these tests that over these four regions the relative differ- 
ences in zero points is less than 0.05 magnitudes for each 
of the four filters. 

4.5. Field galaxy colours 

In this Section we compare our field galaxy colours (or, 
more accurately, extended sources) to those from other 
deep surveys. We compute the median (B — V)ab and 
(V — I)ab colours in half-magnitude bins as a function 
of Bab and Iab sample limiting magnitude. These mea- 
surements are shown as the solid lines in the upper and 
lower panels of Figure El The symbols show field galaxy 
colours measured from the CFDF survey. In both cases 
our agreement with the published values is excellent. 

We also investigate the distribution of galaxies in 
colour-colour space as a function of magnitude; this is il- 
lustrated in Figures ^] ^| and [3U| In these Figures we 
compare our observed galaxy colours to model predic- 
tions for a range of progressively fainter magnitude slices. 
The grey-scales in these Figures are proportional to ob- 
ject number density, computed in 0.1 magnitude inter- 
vals in (B — V)ab and (V — I)ab- Dotted lines show 
the path of late-type (E, SO) galaxies wheras early-type 
galaxies (Sa, Sc and Irregular) are represented by solid 
lines. Galaxy colours were computed at intervals of 0.1 
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Fig. 17. Field galaxy colours from our catalogues. In the 
upper panel, the solid line shows the median (B — V)ab 
galaxy colour computed in half-magnitude bins as a func- 
tion of Bab limiting magnitude; in the lower panel, we 
show the median (V — I)ab colour as a function of 
Iab magnitude. In both cases we compare with colours 
from the CFDF survey (symbols). 

in redshift by applying the CFH12K filter response func- 
tions to model evolving spectral energy distributions gen- 
erated using the "2000" r evision of the GISSEL libraries 
llBruzual fc CharlotL Il99^ . Computations were made as- 
suming a flat A-CDM cosmology. The e-folding times for 
the five tracks are (0.3,1.0,2.0,4.0) Gyrs and "constant" 
for elliptical, SO, Sa, Sc and Irregular galaxy types respec- 
tively. Model galaxy zero-redshift colours (represented by 
the filled d ots) are in good agree ment with the colours 

" ll200ll) Final ly, intergalactic 
Madau (1995) has also been 



reported in iMannucci et al 



absorption as modelled by 
added. 

Figure ^3] shows the bright magnitude slice at 18 < 
Iab < 20. In this Figure, most galaxies lie in a single well- 
defined region. From previous /^-selected galaxy redshift 
surveys such as the CFRS we know that at these magni- 
tudes the median redshift of the faint ga laxy population 
is around z ~ 0.2 i Cramnton et al.lll995h . and this is also 
roughly where the model tracks enter the locus occupied 
by the galaxy population. In addition, we note that the 
model tracks span very well all the observed colours in our 
survey. 

For the intermediate magnitude slice at 20 < Iab < 
22, illustrated in Figure 1191 the (predominantly late- 
type) galaxy population occupies two distinct locii, and in 
colour-colour space there is a comparatively well-defined 
separation between high (z > 0.4) and low (z < 0.4) 
redshift galaxies. For the very faintest magnitude slice, 
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Fig. 18. (B — V)ab vs (V — I)ab colour-colour diagram 
for galaxies with 18 < Iab < 20. The greyscale levels are 
proportional to the numbers of galaxies at each combina- 
tion of (V — I)ab and (B — V)ab colour. The lines shows 
model tracks for five galaxy types (two early types, shown 
by the dotted lines, and three late- types, represented by 
the solid lines). Tracks are computed in the redshift range 
0.0 < z < 0.4 in intervals of 0.1. Solid symbols show the 
location of the z = point for each track. 
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Fig. 19. Colour-colour diagram for galaxies with 20 < 
Iab < 22. For this Figure, model tracks span the redshifts 
range 0.0 < z < 1.0. 



22 < Iab < 24, the sharp distinction has been removed, 
and all parts of the colour-colour diagram are equally 
populated. In general, our extended source colours are in 
agreement with the existing model tracks. 

4.6. Clustering of point-like sources 

In this Section and the following one we investigate the 
clustering properties of point-like and extended sources in 
our catalogues. To do this, we use the well-known pro- 
jected two-point angular correlation function, lo(9). which 
measures the excess of pairs separated by an angle 8,8+69 
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Fig. 20. Colour-colour diagram for galaxies with 22 < 
Iab < 24. Model tracks are computed from 0.0 < z < 3.0. 



with respect to a random distribution. This statistic is use- 
ful for our purposes because it is particularly sensitive to 
any residual variations of the magnitude zero-point across 
our stacked images. 

We first measure the angular correlation function ui(9) 
of the stellar sources. As stars are unclustered, we expect 
that if our magnitude zero-points and detection thresholds 
are uniform over our field then u>(0) should be zero at all 
angular scales. For this test, we consider two magnitude 
limited samples: point-like sources with 18.0 < Iab < 20.0 
and those wi th 20.0 < Iab < 21.5. W e measure uo{9) using 
the standard lLandv fc Szalavl 1 1993h estimator, i.e., 



w(0) = 



DD - 2DR + RR 
RR 



(9) 



with the DD, DR and RR terms referring to the 
number of data-data, data-random and random-random 
galaxy pairs between 9 and 9 + 89. We use logarithmically 
spaced bins with A log(6*) = 0.2. 

Our results are displayed in Figure [2] For both mag- 
nitude cuts, at scales —2.5 < log(#) < 0.0 the measured 
correlation amplitudes are consistent with zero. 

4. 7. Clustering of extended sources 

The angular clustering of field galaxies, parametrised 
by u)(9) provides important information concerning the 
galaxy distribution. Many studies have measured u)(&) for 
progressively larger samples of galaxies rea ching to fainter 
and fa i nter magnitudes (s e e, for example. iBrainerd et alJ 
l|l994lklRoche et alJ l|l993lklEfstathiou et alJjl99lh ). The 
advent of mosaic detector arrays means it is now possible 



Fig. 21. Angular clustering of stellar sources, oj(9) as a 
function of the logarithm of the angular separation in de- 
grees. Filled squares represent objects with 18.0 < Iab < 
20; open circles show sources with 20.0 < Iab < 21.5. For 
clarity, the fainter magnitude slice has been slightly offset 
in log(#) from the brighter one. 



to measure the angular clustering of field galaxies on de- 
gree scales at z ~ 1 with samples reaching ma gnit ude lim- 
its of I ab ~ 24. Both lMcCracken et al.l l)200l|) and lWilsonl 
(2003J) have both measured co(9) using deep catalogues 
constructed using the UH8K camera. In this Section we 
present comparisons between lo{9) measurements in this 
survey with other published works as a further check on 
the reliability of our data reduction procedure. 

In Figure |2U we show the angular correlation func- 
tion uj(9) as a function of angular separation, log(#) for 
19.5 < iAB-med < 23.5 measured from our catalogues. 
We compute lu(9) for -2.8 < log(0) < 0.0 . The solid lines 
show the fitted correlation amplitudes assuming a func- 
tional form u>(9) = A9~ 8 — C parametrisation. The "C" 
term rep resents the in t egral constraint correction (see, for 
example. [Roche et al.1 l|1993ri s |. For our field, we find that 
C = 3.0. Measurements in our catalogue follow the ex- 
pected power-law behaviour at least until angular scales 
where the integral constraint correction becomes impor- 
tant. 

In Figure 12*51 we show the amplitude of lo(9) at 1' (as- 
suming a slope 8 = —0.8) a function Iab sample me- 
dian magnitude for the same magnitude slices presented in 
Figure 12*2*1 The open symbols and s tars show the m easure- 
ments from the CFDF survey and IWilsonl l|2003f) respec- 
tively. The amplitude of oj{9) at 1' steadily decreases with 
limiting magnitude: moreover, all magnitude bins agree 
quite well with the literature measurements. 
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Fig. 22. Jab -selected galaxy clustering. For each magni- 
tude slice, we show the amplitude of the angular correla- 
tion function, u>(9) as a function of the logarithm of the 
angular separation in degrees, log(#). At each magnitude 
slice, we show the fitted correlation amplitude, including 
a correction for the integral constraint term. 



It has recently reported that CFH12K observations 
may be affected by scattered light 7 . This effect origi- 
nates from light reflecting from structures inside the cam- 
era/telescope assembly. It produces an additive correction 
which is not removed by flat-fielding, as the flat-fields 
themselves are not uniformly illuminated. Measurements 
made at CFHT (and described extensively on the ELIXIR 
web pages) show that object magnitudes measured in the 
outer edges of the field of view can be as much 0.1 mag- 
nitudes fainter than in the central region. 

We have looked for possible signatures of this problem 
in our catalogues. At large angular scales, the stellar corre- 
lation function we measure fSection l4.6|) is consistent with 
zero (Figure |2J. However, because the slope of the stel- 
lar number counts is relatively shallow compared to the 
galaxy counts CSection lO|l the radial dependence of stel- 
lar counts will be comparatively weaker (corresponding to 
density changes of around a few percent). For galaxies, the 
steeper count slope means that the radial density depen- 
dence is higher, and so would be easier to detect this with 
a correlation function analysis as outlined above. 

We used a scattered light image supplied by CFHT to 
investigate if the magnitude variation can have a signif- 
icant effect on the measured u>(6). Because of the small 
dither offsets used in our observation strategy we are able 
to separate a single pointing from the four-pointing com- 
bined mosaic, and from this extracted sub-area we 'cor- 



Fig. 23. The amplitude of the angular correlation function 
at arcminute scales. The filled symbols show log A U {V) as 
a function of Iab limiting magnitude as measured in the 
CFH12K-VIRMOS deep field, assuming a power-law form 
for the galaxy angular correlation function uj(9) with slope 
of S = —0.8 (error bars shown represent only the error 
bars from the fitting program and do not indicate the true 
cosmic errors, which are likely to be larger). The open 
circles and stars show measurements from other recent 
wide-field surveys. 

rect' magnitudes in our catalogues based on the radial de- 
pendence derived at CFHT, and then re-apply our selected 
magnitude cuts. For galaxies with 18.5 < Iab < 24.5 we 
find that application of the scattered light map produces a 
difference of around ~ 3.5% in the total numbers of galax- 
ies (below the systematic errors quoted in the previous sec- 
tions). Furthermore, the measured u>(6) values computed 
at each angular separation do not differ significantly before 
and after the application of the scattered light correction 
(the difference is smaller than the Poissonian error bar in 
each bin). From these considerations, we conclude that the 
effects of residual scattered light in our catalogues is not 
detectable, at least to Iab < 24.0. 



5. Summary 

In this paper we have presented BVRI optical data for 
the CFH12K-VIRMOS F02 deep field. Our observations 
reach limiting magnitudes of Bab ~ 26.5, Vab ~ 26.2, 
Rab ^25.9 and Iab ^25.0 (measured from our simula- 
tions). This field consists of a single contiguous area cov- 
ering 1.2 deg 2 . In the magnitude range 18.0 < Iab < 24.0 
there are 90,729 extended sources in our survey. 

We have presented a detailed explanation of the re- 
http://www.cfht.hawaii.edu/Instruments/Elixir/scattered.htidhction and preparation of the object catalogues and dc- 
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scribed how our astrometric and photometric solutions 
were computed. In summary, 68% of catalogue sources 
have an absolute per co-ordinate residual less than |Aa| ~ 
0.38" and |A£| ~ 0.32" in right ascension and declination 
respectively. Expressed in the same sense, our internal 
(filter-to-filter) per co-ordinate astrometric uncertainties 
are |Aa| ~ 0.08" and \A6\ ~ 0.08". Our absolute photo- 
metric calibration, estimated from standard star data, and 
from galaxy number counts, is accurate to ~ 0.1 magni- 
tudes. A series of Monte-Carlo simulations have been used 
to assess the reliability of catalogues in planes defined by 
peak surface brightness, position, and total magnitude. 
These show that our catalogues are reliable, uniform and 
fully complete to Iab < 24.0 . 

Galaxy number counts in BVRI bandpasses have been 
presented. We describe a method of star-galaxy separa- 
tion which correctly accounts for the position-dependent 
seeing present in our final stacked images. Our galaxy 
counts are in good agreement with literature compilations, 
and our stellar counts are well matched by a theoretical 
model of the galaxy. Qualitatively, at each magnitude bin, 
Iab selected galaxy number counts in our survey agree to 
within 10% with measurements made in the CFDF sur- 
vey. Additionally, the slope of the galaxy counts in both 
B and / are consistent with previous determinations from 
other deep surveys. Furthermore, field galaxy colours in 
(V — I) and (B — V) measured in our survey are also 
in agreement (better than 0.1 mag) with measurements 
made in the Canada-France deep fields survey. The loca- 
tion of galaxies in colour-colour space, and its evolution 
with magnitude, is broadly consistent with expectations 
based on simple populations synthesis models. 

We have also measured the angular correlation func- 
tion of both stars and galaxies in our dataset. Stellar cor- 
relations are consistent with zero at all angular scales. 
Clustering of galaxies measured to lAB-med ~ 23 in the 
angular range —3 < log(0) < — 1 follows well the expected 
power law behaviour, and the observed scaling of the am- 
plitude of oj(9) with lAB-med is in agreement with previ- 
ous surveys. We have also investigated the possible effect 
that scattered light could have on the galaxy angular cor- 
relation function, and concluded that to Iab < 24.0 this 
effect does not pose a significant problem to our measure- 
ments. 

Future papers will describe numerous extensions in 
wavelength to the survey field described here. Furthermore 
the catalogues presented in this paper are currently being 
used to select a very large magnitude- limited sample of ob- 
jects for the VIRMOS-VLT deep survey. Together, both 
surveys will provide a unique picture of the Universe at 
intermediate and high redshifts. 
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